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Abstract 

Isoscalar monopole and dipole compressional modes are computed for a vari- 
ety of closed-shell nuclei in a relativistic random-phase approximation to three 
different parametrizations of the Walecka model with scalar self- interactions. 
Particular emphasis is placed on the role of self-consistency which by itself, 
and with little else, guarantees the decoupling of the spurious isoscalar-dipole 
strength from the physical response and the conservation of the vector current. 
A powerful new relation is introduced to quantify the violation of the vector 
current in terms of various ground-state form-factors. For the isoscalar-dipole 
mode two distinct regions are clearly identified: (i) a high-energy component 
that is sensitive to the size of the nucleus and scales with the compressibility 
of the model and (ii) a low-energy component that is insensitivity to the nu- 
clear compressibility. A fairly good description of both compressional modes 
is obtained by using a "soft" parametrization having a compression modulus 
of if = 224 MeV. 
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I. INTRODUCTION 



The study of nuclear compressional modes, while interesting in its own right, is motivated 
by our desire to understand the equation of state of hadronic matter, especially in relation 
to its compression modulus. In turn, an accurate determination of the compression modulus 
places important constraints on theoretical models of nuclear structure, heavy-ion collisions, 
neutron stars, and supernovae explosions. 

While it remains true that measuring the energy of the nuclear compressional modes pro- 
vides the most accurate determination of the compression modulus, significant advances in 
astronomical observations and terrestrial experiments are providing important complimen- 
tary information. For example, explaining the time structure of the neutrino burst emitted 
from supernova SN1987A seems to require a relatively soft equation of state as input in the 
simulations of core-collapsed supernova Further, the recently inferred narrow mass 

distribution of neutron stars |J poses stringent constraints on the nuclear equation of state. 
At the same time, a number of improved radii-measurements of radio-quite, isolated neutron 
stars — such as RX J185635-3754 — will contribute significantly to our understanding of the 
high-density component of the equation of state ||]. Finally, measurements of the elliptical 
flow in relativistic heavy-ion reactions seem to have established the utility of this observable 
as a probe of the stiffness of the equation of state H. 

Also significant is the strong correlation between seemingly unrelated experiments. In- 
deed, the radius of a neutron star is predicted to be strongly correlated to the neutron skin 
of a heavy nucleus ||[7| . Thus, the upcoming measurement of the neutron radius of 208 Pb at 
the Jefferson Laboratory IIU should place important limits on the radii of neutron-stars. 

Although measurements of the giant monopole resonance [10,11] and the isoscalar giant 
dipole resonance |l2|-p^] have existed for some time, the field has seen a revitalization due 



to new and improved measurements of both compressional modes fli5|-|r7j| . The field has 
also seen significant advances in the theoretical domain. Indeed, calculations of nuclear 
compressional modes using Hartree-Fock (HF) plus RPA approaches with state-of-the-art 
Skyrme interactions are now possible []T8| , [T9| . Relativistic RPA models have also enjoyed a 
great deal of success, especially now that scalar self-interactions have been incorporated into 
the calculation of the response p0|-[23|| . At the same time the philosophy behind the theo- 
retical extraction of the nuclear compressibility has evolved considerably. Earlier attempts 
depended heavily on semi-empirical formulas that related the compressibility to the ener- 
gies of the compressional modes [24]. The field now demands stricter standards: the model, 
without any recourse to semi-empirical mass formulas, must predict both the compressibility 
of nuclear matter as well as the energy of the compressional modes. 

In this publication state-of-the-art calculations of the isoscalar giant- monopole resonance 
(GMR) and the isoscalar giant-dipole resonance (ISGDR) are reported for a variety of closed- 
shell nuclei. This paper represents an expanded version of a short article published recently 
that focused exclusively on 208 Pb [p2] . The model adopted in this work is based on a 



relativistic random-phase-approximation (RPA) to three different parameterizations of the 
Walecka model with scalar self-interactions. A nonspectral approach that treats discrete 
and continuum excitations on equal footing is implemented. As a result, the conservation 
of the vector current is strictly maintained throughout the calculation. Moreover, for the 
calculation of the RPA response we employ a residual particle-hole interaction consistent 
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with the particle-particle interaction used to generate the mean-field ground state. In this 
way the spurious isoscalar-dipole strength, associated with the uniform translation of the 
center-of-mass, gets shifted to zero excitation energy and is cleanly separated from the 
physical response. 

Having established the theoretical underpinning of our calculation, it is now useful to 
contrast it against alternative self-consistent implementations. In a recent article by Shlomo 
and Sanzhur [25 1, it is suggested that actual implementations of the RPA, in spite of claim- 
ing otherwise, are not fully self-consistent. It is pointed out that these calculations often 
resort to a variety of approximations such as: (i) neglecting the two-body Coulomb and 
spin-orbit terms in the residual particle-hole interaction, (ii) approximating the momentum- 
dependent parts in the particle-hole interaction, (iii) limiting the particle-hole space in a 
discretized calculation by a cut-off energy -Eyf x , and (iv) introducing a smearing parameter, 
such as a Lorentzian width. Each of these approximations is now briefly addressed. In the 
relativistic formalism employed here neither the two-body Coulomb nor the spin-orbit inter- 
action are neglected. Rather, the residual particle-hole interaction includes the (isoscalar) 
contribution from the photon as well as spin-orbit effects that are incorporated — to all 
orders — by merely maintaining the relativistic structure of the interaction. Moreover, the 
residual particle-hole interaction is momentum independent because one preserves intact 
its full Lorentz structure; no momentum-dependence is generated through a nonrelativistic 
reduction of the interaction. Further, the non-spectral approach employed here avoids any 
reliance on artificial cutoffs and truncations. Finally, while a Lorentzian width is included to 
compute the properties of discrete excitations, it is done so by ensuring that the physically 
relevant quantities, the excitation energy and the inelastic form-factor, remain invariant 
under a change in width. 

The paper has been organized as follows. Section |TJ describes the relativistic mean-field 
plus RPA formalism in great detail placing special emphasis on the role of self-consistency. 
Section [nj illustrates the importance of self-consistency for the conservation of the vec- 
tor current and for the decoupling of spurious strength from the physical isoscalar-dipole 
response. Here a powerful novel relation is introduced to quantify the violation of the vec- 
tor current in terms of various known ground-state form-factors. Results are displayed in 
Sec. [IV], while a summary and conclusions are presented in Sec. M. 



II. FORMALISM 

In this section a detailed description of the mean-field plus RPA formalism employed 
to compute the distribution of strength for both compressional modes is presented. This 
formalism, with the exception of its implementation in the case of scalar self-interactions, 
has now been available for almost fifteen years [p6j -28[. However, important lessons keep 



being ignored [pD[ , just to be soon rediscovered pif . Thus, we feel compelled to present, for 
what we hope is the last time, a thorough discussion of the relativistic RPA formalism. 

The first step in calculating a relativistic RPA response is the computation of the mean- 
field ground state in a self-consistent approximation. Once self-consistency is achieved, 
three important pieces of information become available: (i) the single-particle energies of 
the occupied orbitals, (ii) their single-particle wave functions, and (iii) the self-consistent 
mean-field potential. This mean-field potential, without any modification, must then be 
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used to generate the nucleoli propagator; in this way the conservation of the vector current 
is guaranteed to be maintained. The nucleon propagator is computed nonspectrally to 
avoid any dependence on the artificial cutoff's and truncations that plague most spectral 
approaches. Moreover, through a nonspectral approach one gives equal treatment to both 
bound and continuum orbitals. 

Having generated the occupied single-particle spectrum and the nucleon propagator, the 
computation of the lowest-order (Hartree) polarization is reduced to the evaluation of various 
matrix elements of the relevant transition operator. To compute the RPA response one needs 
to go beyond the single-particle response. The RPA builds coherence among the many 
allowed particle-hole excitations by iterating the lowest-order polarization to all orders via 
the residual particle-hole interaction. Yet special care must be taken in adopting a residual 
particle-hole interaction consistent with the particle-particle interaction used to generate 
the mean-field ground state. Only then can one ensure that the spurious component of the 
isoscalar-dipole response will get shifted to zero excitation energy P^| , |SD[| . As the polarization 
tensor is a fundamental many-body operator, it can be computed systematically using well- 
known many-body techniques pi| . Having computed the polarization tensor, the nuclear 
response is extracted by simply taking its imaginary part. The following sections provide a 
detailed account on the implementation of these ideas. 



A. The Lagrangian Density 



The starting point for the calculation of the nuclear response is a Lagrangian density 
having an isodoublet nucleon field (if)) interacting via the exchange of two isoscalar mesons, 
the scalar sigma (0) and the vector omega (V M ), one isovector meson, the rho and the 
photon (A**) 
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The pseudoscalar pion is not included as it does not contribute at the 
mean-field level. In addition to meson-nucleon interactions the Lagrangian density includes 
scalar self-interactions. These are responsible for reducing the nuclear compressibility from 
the unrealistically large value of K = 545 MeV, obtained in the original linear model of 
Walecka j34||, all the way down to the acceptable value of K = 22A MeV. Thus, without the 
inclusion of scalar self-interactions a realistic calculation of the compressional modes is not 
feasible. The Lagrangian density for the model is thus given by 

Ant = $[gd -g v y-^T-f>- ~(1 + t 3 )A}^ - U(<f>) ; U{<f>) = ^ + iA0 4 , (1) 

were use of the "slash" notation, Y= 7 M ^> has been made. The various model parameters 
have been listed in Table ID. 



B. The nucleon propagator 

The mean-field propagator contains information about the interaction of the propagating 
nucleon with the average potential generated by the nuclear medium. However, even in a 
Fermi-gas description, where all interactions are neglected, the nucleon propagator would 
still differ from its free-space value because of the presence of a filled Fermi sea. Indeed, 
the analytic structure of the free-nucleon propagator at finite density is different from its 
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free-space value (see Fig. |l|). This suggests the following decomposition of the nucleon 
propagator |32| : 



G(x,y) 



du 
— < 
2tt 



-iui(x°— y°) 



G(x,y; u) 



G(x, y; oj) = G F (x, y; uj) + G D (x, y; uj) . 



(2a) 
(2b) 



The Feynman part of the propagator, Gp, admits a spectral decomposition in terms of the 
mean-field solutions to the Dirac equation. That is, 



G F (x,y;u) = 



U n (x)U n (y) + V n {x)V n (y) 



UJ - E n +) + 17] UJ + En 



(") 



IT] 



(3) 



where U n and V n are the positive- and negative-energy solutions to the Dirac equation, and 
the sum is over all states in the spectrum. The analytic structure of Gp is identical to that of 
the conventional Feynman propagator [35]. The density- dependent part of the propagator, 
Gd-i corrects Gf for the presence of a filled Fermi sea. This correction occurs even in a 
noninteracting system and is due to the Pauli exclusion principle. Formally, one effects this 
correction by shifting the position of the pole of every occupied state from below to above 
the real axis (see Fig. [I]) 



G D (x,y;^)= £ C/ n (x)C/ n (y) 



n<F 



IT] 



En +) + if]. 



2mJ2 S(u - EM)u n (x)U n (y) 



(4) 



n<F 



Note that the sum over n is now restricted to only those positive-energy states below the 
Fermi energy. In a mean-field approximation these states satisfy a Dirac equation of the 
form: 



£i +) 7° + *7 ■ V - M - Emf(x)] C/ n (x) = 
where the mean-field potential is given by 

Smf(^) = S s (x) +7°S (x) . 



(5) 



(6) 



The quantities Ss and So denote the scalar and vector potentials that have been generated 
self-consistently at the mean-field level. Since this work is limited to the response of closed- 
shell nuclei, it is assumed that the mean-field potential has been generated by a spherically- 
symmetric, spin-saturated ground state. 

Although the above spectral decomposition of the nucleon propagator will become im- 
portant in understanding the spectral content of the nuclear response, in practice it suffers 
from a reliance on artificial cutoffs and truncations. An efficient scheme that avoids such 
a dependence is the nonspectral approach. A nonspectral approach has the added advan- 
tage that both positive- and negative-energy continuua are treated exactly. As a result, the 
contributions from the negative-energy states to the response are included automatically. 
This is important to maintain fundamental physical principles, as the positive-energy states 
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by themselves are not complete. To obtain the nucleon propagator in nonspectral form one 
must solve the following inhomogeneous Dirac equation: 



u; 7 ° + *7 • V - M - £mf(z)1 G f (x, y; u) = 5(x - y) . 



(7) 



Here w is taken to be a complex variable and the mean-field potential is identical to the 
one used to generate the nuclear ground state. Taking advantage of the spherical symmetry 
of the potential, one may decompose the Feynman propagator in terms of spin-spherical 
harmonics 



G F (x,y;u) = — X] ( 

Km \ 



/ gi 1 (x,y;u)('k\ + Km)(+K,rn\y) —igi 2 (x,y;w){k\ + Km)(—Km\y) 
. ig-tiix, y, u) (xj -nm) (+nm\y) g$ 2 (x, y; u) (xj - nm) (-nm\y) 



which are defined as 



(x|«m) = Yl ( lm lA m s\Hj m ) Y lrni(x)Xh 



raim a 



3 = \ K \ 



- and I = 
2 



+ K if K > , 

— K — 1 if K < . 



(8) 

(9) 
(10) 



The above decomposition enables one to rewrite the Dirac equation as a set of first-order, 
coupled, ordinary differential equations of the form 



uj*-M* 



V dx x 



d k* 
dx x 
-lu*-M* 







) U21 


922) 



8(x - y) 



(11) 



where we have defined 



EJx) and M* = M + E a (x) . 



(12) 



It is important to underscore that the mean-field potentials used to compute the nucleon 
propagator must be identical to those used to generate the mean-field ground state if the 
conservation of the vector current is to be maintained. 



C. The nuclear polarization 

To illustrate the many-body techniques employed in the manuscript, we define a general 
polarization insertion as the time-ordered product of two arbitrary nucleon currents: 

*o> , (13) 

where ^ denotes the exact nuclear ground state and J a (x) is a one-body current operator 
of the form 

J a (x) = i)(x)Y a i)(x) . (14) 



iU a ^(x,y) = {$ \T J a (x)JP(y) 
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Note that the "big" gamma matrices have been defined so that the one-body current operator 



be hermitian [35]. That is, 

r a = {i, rf, y, 7 V, (t^} with r c 



7 ° r Qt 7° 



■pa 



(15) 



In a mean-field approximation to the nuclear ground state, such as the one employed here 
and in most of the other relativistic calculations to date, the polarization insertion may be 
written exclusively in terms of the nucleon mean-field propagator 



Il a P(x,y) = Ti\r a G{x,y)T p G{ 



(16) 



The earlier decomposition of the nucleon propagator into Feynman and density- 
dependent contributions [Eq. (|2b|) ] suggests an equivalent decomposition for the polarization 
insertion 



J-oo 2n 

IT* (x, y; lj) = Uf (x, y; u>) + 11? (x, y; u) . 



(17a) 
(17b) 



The Feynman part of the polarization, IT^ 3 , is independent of Gd and describes the po- 
larization of the vacuum. This piece, which diverges and needs to be renormalized, has 
been incorporated in our earlier calculations of the longitudinal response in the quasifree 
region p6|| . However, it has been included only in a local-density approximation. To our 
knowledge an exact finite-nucleus calculation of vacuum polarization has yet to be per- 
formed. While a local-density approximation is accurate in the quasifree region where many 
angular-momentum channels contribute, it has proven inadequate for the description of dis- 



crete nuclear excitations 37 ] . In particular, the spurious isoscalar dipole strength associated 
with the uniform translation of the center of mass does not get shifted all the way down to 
zero excitation energy. More relevant, the role of vacuum polarization in effective hadronic 
theories is currently being revisited. Effective Field Theories now suggest that the largely 
unknown physics associated with the short- distance dynamics may be effectively simulated 
by the use of various local operators []3"5|-fi0|. It is for these reasons that vacuum polarization 
will be ignored henceforth. Note, however, that it is still possible to ignore vacuum effects 
and end up with a completely consistent model of the nuclear response [|3"D,3"2||. 

In contrast to the Feynman part of the polarization, the density-dependent part is finite 
and can be computed exactly in the finite system |26|-p8|. It is given by 



n^f(x, y; u) = n^(x, y; u) + IT^(x, y; 



where 



ri^(x,y;u,) = ]T f/ n (x)r Q G F (x,y;+a; + E^)v^U n {y) 



n<F 



(18) 

(19a) 
(19b) 



n<F 



Note that the Pauli-blocking of particle-hole excitations, a term usually denoted by 11^, 
has already been incorporated in the above two terms. The density-dependent part of the 
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polarization includes the excitation of particle-hole pairs plus the mixing between positive- 
and negative-energy states; this last term is sometimes referred to as the Pauli blocking 
of NN excitations. The spectral content of II^ is easily revealed by using the spectral 
decomposition of the Feynman propagator [see Eq. (§)]. For example, the Feynman-density 
component of the polarization, IT^,, may be written as 

I7 n (x)r^ m (x)F m (y)r^ n (y) | F n (x)r"K re (x)F m (y)r^ n (y) " 
uj - (MP - eL +) ) +i V u+ (E^ ] + E { „ +) ) -irj 

(20) 

The first term in the sum represents the excitation of a particle-hole pair. The excitation 
becomes real, namely both particles go on-shell, when the energy transfer to the nucleus 
becomes identical to the pair-excitation energy uj = E^' — E^K The second term in the sum 
has no nonrelativistic counterpart; it represents the mixing between positive- and negative- 
energy states. Although the contribution from vacuum polarization has been neglected, 
the inclusion of this mixing is of utmost importance for maintaining current conservation. 
Moreover, it is also essential for the removal of all spurious strength from the excitation of 
the isoscalar dipole mode. The inclusion of the negative-energy sector in the calculation of 
the response underscores the basic fact that the positive-energy sector of the spectrum, by 
itself, is not complete. 



n^j(x,y;w)= Yl 

m,n<F 



D. The RPA equations 

The polarization tensor describes modifications to the propagation of various mesons 
(such as the a, uj, p, . . .) as they move through the nuclear environment. In addition, 
the polarization tensor contains all information on the excitation spectrum of the nucleus. 
Indeed, the polarization insertion is an analytic function of the frequency uj, except for the 
presence of simple poles located at the excitation energies of the system. The residue at the 
pole is simply related to the inelastic form-factor [[31]. 

The singularity structure of the lowest-order polarization tensor is easily inferred from the 
mean-field spectrum: the nuclear excitation energies (poles) appear at energies given by the 
difference between the single-particle energies of a nucleon above the Fermi level (particle) 
and one below (hole). In this approximation the residual interaction between the particle 
and the hole is neglected. However, the consistent response of the mean-field ground state 
demands that the residual interaction between the particle and the hole be incorporated [|30 



This may be implemented by solving Dyson's equation for the polarization insertion in a 
random-phase approximation. In RPA the lowest-order polarization is iterated to all orders 
via the residual particle-hole interaction. Because the iteration is to all orders, the singularity 
structure of the propagator, and thus the location of the poles, is modified relative to the 
lowest-order predictions. Dyson's equation for the RPA polarization is given by: 

d 3 k d 3 k' 
(2tt) 3 (2tt) 

where V\ a (k, k'; uj) is the residual interaction to be discussed below and Il^f (q, q'; u) is the 
Fourier transform of the lowest-order polarization. That is, 



r d 3 k cfik' 

nf PA (q, q'; uj) = IT^(q, q'; uj) + y _ ™ I#( q , k; uj)V Xa (k, k'; uj)U^ PA (k f , q'; uj) , (21) 



S 



n£f (q, q'; u>) = J d 3 x d 3 y e^-W-y) n^(x, y; w ) . (22) 

At this point it is convenient to depart from the general formalism adopted until now and 
restrict the discussion to the case of interest: the isoscalar compressional modes. Hence, the 
only component of the residual interaction that must be retained is the one mediated by 
the exchange of the sigma and omega mesons, and the (isoscalar component of the) photon. 
Moreover we employ the simplest operator, the timelike component of the vector current 

p(q) = f d 3 a;e iq - x ^(x)7V(x) , (23) 

that can couple to these natural-parity excitations. 

The computational demands imposed on a calculation of the RPA response for a nucleus 
as large as 208 Pb can be formidable indeed. Powerful symmetries that are present in infinite 
nuclear matter, such as translational invariance, are broken in the finite system. As a result, 
the RPA equations that were algebraic in the infinite system become integral equations 
in the finite nucleus. Moreover, modes of excitation that were uncoupled before, such as 
longitudinal and transverse modes, become coupled now. In this way the RPA equations, 
because of the ubiquitous scalar-longitudinal mixing, become a complicated 5x5 set of 
coupled integral equations. Correspondingly, the residual particle-hole interaction, also a 
5x5 kernel, may be written as 

v (V \c'- f A - ^s A (k,k';a;) \ , . 

v ^ k ^-{ ^(k,k»J ' 

where the vector propagator is given by 



D aP (k, k'; W ) = (2vr) 3 5(k - k') (- 9aP + ) D(k, u) ; D(k, u) = . (25) 

mi ) u z — k z — mi 



Note that because vector self-interactions have not yet been included in the present version 
of the model, the vector propagator remains local (in momentum space) and maintains its 
simple Yukawa form. In contrast, scalar self- interactions modify the propagator relative to 
its simple free-space form. Hence, the scalar propagator now satisfies a nontrivial Klein- 



Gordan equation of the form: [20,21 



^ + V 2 -mi-U"(<j ) ) A(x,y; W ) = *(x-y) . (26) 



III. FUNDAMENTAL SYMMETRIES 

In the following two sections we discuss important symmetries related to the conservation 
of the vector current and to the elimination of the spurious isoscalar-dipole strength from 
the physical response. We are adamant about the preservation of these two fundamental 
symmetries of nature as we regard the predictions of theoretical formulations that violate 
them as ambiguous at best. For example, in a framework that violates the conservation of 
the vector current should one calculate the longitudinal response of the nuclear ground state 
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by using the timelike component or the longitudinal one? Likewise, the predicted energy 
and distribution of isoscalar-dipole strength in a model that retains even a small fraction of 
spurious strength will bear little resemblance to reality. It is only through consistency, the 
recurring theme of this paper, that one can enforce these important dynamical demands. 
How is that consistency plays such an important role in preserving these fundamental sym- 
metries, will now be discussed. 



A. Conservation of the vector current 

We start by discussing the conservation of the vector current. Current conservation 
demands that the timelike component of vector current be related to the longitudinal com- 
ponent. This impacts greatly on the results; it forces the nuclear polarization with one 
Lorentz vector index to be transverse to the four-momentum transfer, irrespective of the 
Lorentz character of the other vertex. That is, 

g Ai n^(q,q» = with g"=(w,q). (27) 



So how is current conservation realized in our model? As indicated in Eq. (|18|) the density 
dependent part of the polarization tensor consists of two terms: Tlp D and ffj^. Does 
each term separately satisfy current conservation or does the conservation of the current 
depend on a sensitive cancellation between them? To address this question we introduce 
the longitudinal (with respect to q) component of the vector current. We start with the 
Feynman-density piece: 

gir^(q,q';a,)= / d 3 xd 3 y £ Z7 n (x)( 7 ■ q)e^ x G F (x, y;uj + f)e^ T^U n (y) . (28) 

J n<F 

To make contact with the timelike component of the polarization we turn the momentum 
transfer q into a gradient operator [(7 • q)e _:(C| X = (7 • zV)e _iq x ] and integrate by parts. 
In this way the gradient operator acts now on both the bound-state nucleon spinor and 
the nucleon propagator. It is then the difference between their respective Dirac equations 
[Eqs. (^) and (0)] that dictates how severe the violation of current conservation becomes. 
We obtain 

g M n^(q,q» = /(q-q') = / d^xe'^^ £ F n (x)r% n (x) , (29) 

J n<F 

where pr (q) represents a ground-state form-factor. This is a new and important result. First, 
such a simple relation would have been impossible to obtain had the mean-field potential for 
the nucleon propagator been any different than the corresponding one for the bound-state 
wave function. This is one of the many manifestations of consistency in the formalism. 
Second, because in spherical nuclei all form-factors are real [32], the imaginary part of IT^,, 
by itself, satisfies current conservation. However, this is not true for the real part. Indeed, 
the violation to the real part of the polarization is regulated by the various ground-state 
form-factors. This result may be used as a stringent test on the numerics. For instance, if 
one lets —>■ 7 and sets q = q' in Eq. (p9[) , the violation becomes identical to the mass 
number of the nucleus. That is, 
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g M nf D (q,q;a;) = / ' d?x £ [/ n (x) 7 ^,(x) = / d 3 a;pB(x) = A . (30) 

^ n<F J 

In Fig. ^ we display the cumulative violation of the vector current as a function of the 
angular-momentum channel J w . Note that the plot also includes the corresponding violation 
in the Density- Feynman part of the nuclear polarization which is given by 

3A(q,q'i W ) = V(q-q')- (31) 
In this way current conservation is properly restored: 



^Il£f (q, q'; u) = q, [Uf D (q, q'; u) + II&(q, q'| w)J = . (32) 

Note that current conservation is maintained for each individual J^-channel. Figure [| vali- 
dates this statement by displaying the timelike component of the polarization alongside the 
longitude component (3 = q) for the isoscalar-dipole state in 40 Ca. These results emerge 
from two powerful demands. First, the interaction driving the nucleon propagator must be 
identical to the one generating the mean-field ground state. Second, the negative-energy 
part of the spectrum must be kept, otherwise the nucleon propagator fails to become the 
Green's function for the relevant Dirac problem. One of the great virtues of the nonspectral 
approach is that the negative-energy states are included automatically. 

So far our discussion of current conservation has been limited to the lowest-order po- 
larization. Nevertheless, the conservation of the vector current at the RPA level places no 
additional demands on the formalism. Indeed, it relies exclusively on the conservation of the 
vector current at the Hartree level and it is independent of the nature of the residual inter- 
action. This result may be derived from the structure of Dyson's equation for the nuclear 
polarization. Using Eqs. ([H]) and (|32"D we obtain 



^ n RPA(q 5 q'; = g M n£f (q, q'; <j) 

d 3 k d 3 k' 



+ 



(2tt) 3 (2tt) 3 



q.U^iq, k; u) V XrT (k, k'; u)U^ PA (k' , q'; to) = . (33) 



We close this section with a brief comment. As the conservation of the vector current is exact 
in our formalism, we are entitled to a minor simplification: the longitudinal component of 
the current can be systematically eliminated in favor of the timelike component. Thus, the 
RPA equations may be reduced from a 5 x 5 to a 4 x 4 set of integral equations by simply 
adopting a modified longitudinal propagator of the form: 

D (k,Lu) = (^\ D(k,u) ; kl = (oo 2 -k 2 ) . (34) 

Note that the gauge component of the vector propagator [the k a kp term in Eq. (p5|) 1 has 
been eliminated from any further discussion because the vector mesons do indeed couple to 
a conserved vector current. 
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B. Spurious strength in the isoscalar-dipole response 



While we have argued earlier that the conservation of the vector current at the RPA level 
is maintained irrespective of the nature of the residual particle-hole interaction, a consistent 
residual interaction becomes of utmost importance in the elimination of the spurious strength 
from the isoscalar-dipole response. This result, first demonstrated by Thouless for the 
nonrelativistic case [M and later extended by Dawson and Furnstahl to the relativistic 



domain |3(| , reinforces the importance of consistency in the formalism. As in the case of the 
conservation of the vector current, the decoupling of the spurious component of the isoscalar- 
dipole response depends on the consistency between the residual particle-hole interaction 
and the particle-particle interaction driving the mean-field ground state. Figure |], where 
the distribution of isoscalar-dipole strength in 16 is displayed, elucidates this point in a 
particularly clear fashion. The lowest-order Hartree response (dashed line) concentrates 
most of the isoscalar-dipole strength in a single fragment located around lo — 16 MeV of 
excitation energy. This is the region where many single-particle transitions from the p- 
shell to the sd-shell occur. Yet most of this strength is spurious, as evinced by the large 
amount being shifted to zero excitation energy in the RPA response (solid line). What 
remains is a relatively small fragment centered around u = 10 MeV of excitation energy; 
we identify this fragment as the first physical isoscalar-dipole state in 16 O. We have also 
included in Fig. | an RPA calculation (dot-dashed line) with a slightly "tampered" residual 
interaction, namely, one that neglects the contribution from the isoscalar component of the 
photon. Although much weaker than its purely isoscalar (sigma and omega) counterparts, 
the photon contribution remains indispensable at low-excitation energies. Indeed, without 
it the spurious center-of-mass state fails to move all the way down to zero excitation energy. 

A similar calculation for the linear L2-set is displayed in Fig. [|. This time, however, 
the width has been reduced considerably (from rj — 1 MeV to i] = 0.05 MeV) so that the 
various discrete single-particle excitations (dashed line) may be resolved. For example, the 
two small fragments in the 10-12 MeV region (dashed line) represent the proton and neutron 
IP 1 / 2 — » 2S l l 2 single-particle excitations respectively (see Table |H[). Moreover, by reducing 
the width one removes any contamination from the spurious state into the first physical 
excitation (solid line). This is essential for a reliable extraction of the inelastic form-factor, 
which is proportional to the area under the peak: 



Fl(q) = \im— S L (q,uj)du. (35) 



Here uj n represents the (discrete) excitation energy. In Fig. |6] we show the isoscalar dipole 
form-factor extracted from the longitudinal response. As we compare with actual experi- 



mental data [fy]], the single-nucleon form-factor has been folded into the calculation. The 
Hartree form-factor is the Fourier transform of the IP 1 / 2 — > 2S* 1 / 2 single-particle transition 
density. As such, it displays a very deep minimum due to the presence of a node in the 2S 1 ^ 2 
wave function. Clearly, even a small amount of configuration mixing will fill in this minimum. 
Indeed, not only does the RPA form-factor (solid line) shows no evidence of a minimum, but 
it actually peaks very close to the Hartree minimum. Further, if the separation between the 
spurious state and the physical states is complete, then the momentum-transfer dependence 
of the isoscalar dipole form-factor should display an octupole (J = 3) behavior rather than 
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that of a dipole |30| , |37| , ^2| . It may be seen in Fig. || that the g-dependence of the physical 
form-factor is indeed (practically) identical to that of the octupole form-factor. 



IV. RESULTS 

Having established the theoretical framework for the calculations of the response, we 
now proceed to display our results for the distribution of isoscalar monopole and isoscalar 
dipole strength on a variety of closed-shell nuclei. As both monopole and dipole states can 
be excited through the timelike component of the vector current, we limit our discussion to 
the longitudinal response: 

5 L (q,w) = ^|(^„|p(q)|^o)| 2 <K^-^) 

n 

= -ij m n 00 (q,q,u;), (36) 

7T 

where p(q) is the Fourier transform of the isoscalar vector density, \I/o is the exact nuclear 
ground state, and ^f n is an excited state with excitation energy uj n . 



A. Isoscalar Giant Monopole Resonance 

The isoscalar giant monopole resonance is the quintessential compressional mode. Re- 
garded as the "breathing mode" of the nucleus, this excitation holds a special place in nuclear 
physics as it provides, perhaps more than any other measurable observable, the most direct 
determination of the compressibility of nuclear matter. In Fig. |7] we display the distribution 
of isoscalar monopole strength for various closed-shell nuclei as predicted by three relativis- 
ts mean-field models. These predictions are also compiled in Table |TTT| . The three models 
have been defined in Ref. [33| as L2, NLB, and NLC and have been constrained to repro- 



duce several bulk properties of nuclear matter at saturation as well as the root-mean-square 
charge radius of 40 Ca; the last two models include self-interactions among the scalar field. 
The model parameters have been listed in Table |. First discovered in a-scattering exper- 
iments from 208 Pb |10| and recently measured with better precision, the peak of the GMR 



has been reported to be located at an excitation energy of E— 14.2±0.1 WJ\. As reported 
in a recent publication [p^j , we found reasonable agreement between experiment and our 
theoretical calculations using set NLC. The other two sets, with compression moduli larger 
than K = 420 MeV, predict the location of the GMR at too large an excitation energy. This 
behavior continues all through the periodic table. Indeed, for medium-size nuclei, such as 
16 O and 40 Ca, it becomes difficult to even identify a genuine GMR with parameter sets L2 
and NLB. In contrast, the identification of the GMR with parameter set NLC is unambigu- 
ous for all nuclei and its prediction for the location of the GMR in 90 Zr is in good agreement 
with experiment. Finally, acceptable agreement has been found with empirical formulas that 
suggest that the position of the GMR should scale as the square root of the compressibility. 
For example, peak energies for this mode have been computed in the ratio of 1:1.38:1.53 for 
208 Pb and 1:1.43:1.57 for 90 Zr, while the square root of the nuclear-matter compressibilities 
are in the ratio of 1:1.37:1.56. These results suggest that models of nuclear structure having 
compression moduli well above K^200 MeV are likely to be in conflict with experiment. 
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B. Isoscalar Giant Dipole Resonance 



The special role played by the isoscalar giant dipole resonance in constraining theoretical 
models of the nuclear response has been discussed extensively in previous sections. We trust 
that our results have convinced the reader that the approach is sound and that the spurious 
contamination has been efficiently removed from the physical excitations. Hence, in the 
remainder of this section we focus on the nuclear and model dependence of the ISGDR. 
Moreover, we also discuss the substantial amount of isoscalar dipole strength predicted to 
exist at low energy and already observed experimentally. 

The distribution of isoscalar-dipole strength in 90 Zr and 208 Pb is displayed in Fig. ^]for the 
three relativistic models. Note that no attempt has been made to identify a resonant peak 
for the case of 16 and 40 Ca as the strength becomes too fragmented. As remarked earlier, 
the model with the softest equation of state (NLC) provides the best description of the 
experimental data P2| . Thus, as hoped, the high-energy component of the isoscalar-dipole 
response provides an independent determination of the compression modulus of nuclear 
matter. Moreover, it constraints, more than any other observable, theoretical models of the 
nuclear response. Even so, we should note that the most accurate of the models (NLC) still 
overpredicts by almost 5 MeV the energy of the isoscalar dipole mode in 90 Zr. 

In contrast to the high-energy component of the isoscalar-dipole response, the low-energy 
component is independent of the compression modulus of nuclear matter (see Fig. ||). Indeed, 
the lowest-energy fragment in 208 Pb is located at an excitation energy of about 8 MeV - 
irrespective of the parameters of the model. That is, relativistic models having compression 
moduli ranging from 220 MeV all the way up to 550 MeV predict a similar distribution 
of low-energy isoscalar-dipole strength in 208 Pb. This behavior continues all throughout 
the periodic table. While the extraction of a sole RPA state, and thus of an associated 
form-factor, is difficult in the case of heavy nuclei, some interesting features emerge from 
the study of the momentum-transfer dependence of the distribution of strength. Figure [L0 
displays such a dependence for 208 Pb. It shows that the large amount of spurious strength 
observed at low-momentum transfer (q = 45 MeV) in the Hartree response gets shifted to 
zero excitation energy (not shown in this figure) leaving a barely visible physical fragment 
at around 7-8 MeV. Moreover, the evolution of RPA strength with momentum transfer seem 
to follow the trends displayed by the inelastic form factor of 16 (see Fig. ^). It has been 
proposed in Ref . |23| , from an analysis of the velocity fields, that the low-energy component 



of the isoscalar dipole mode is determined by surface effects. 



V. CONCLUSIONS 

The distribution of isoscalar monopole and isoscalar dipole strength has been computed in 
a relativistic random-phase-approximation to the Walecka model using various parametriza- 
tions that incorporate scalar self-interactions. While all of these models provide an equally 
good description of the properties of nuclear matter at saturation, their predictions for the 
nuclear compressibility differ by more than a factor of two. Predictions for the energy 
of various surface modes in medium-mass nuclei within a self-consistent random-phase- 
approximation to the Walecka model have existed for over a decade. However, attempts 
at calculating compressional modes in the original model, with a compressibility of K = 545 
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MeV, were doomed to failure. Recently, however, scalar self-interactions, so instrumental in 
softening the equation of state, were incorporated into the calculation of the response. Un- 
fortunately, as lessons were being learned, others were being forgotten. Chief among these 
was the important role of the negative-energy states in the formalism. 

In this paper the relativistic RPA formalism with scalar self-interactions has been re- 
viewed in great detail. A nonspectral approach has been implemented that automatically 
includes both positive- and negative-energy continuua without any reliance on artificial cut- 
offs and truncations. Special emphasis was placed on the role of self-consistency which 
demands that the same interaction used to generate the mean-field ground state be used to: 
(i) compute the nucleon propagator and (ii) the RPA response. Enforcing (i) guarantees the 
conservation of the vector current, while enforcing (ii) successfully decouples the spurious 
isoscalar-dipole strength from the physical response. A novel relation that quantifies the vio- 
lation of the vector current exclusively in terms of ground-state form-factors was introduced. 
This relation may be used as a stringent test on the numerics. 

Predictions for the isoscalar giant-monopole resonance in the NLC model, with a nuclear 
compressibility of K = 224 MeV, were in good agreement with experiment and also with 
semi-empirical formulas that suggest that the position of the GMR should scale as the square 
root of the nuclear compressibility. For the isoscalar-dipole mode the best description of the 
data was still obtained with the NLC set, but here the discrepancies were larger than in 
the monopole case. In particular, theoretical calculation overestimate the position of the 
ISGDR in 90 Zr by almost 5 MeV. In addition to the high-energy component of the ISGDR, 
a low-energy component that is insensitive to the compressibility of the model was clearly 
identified in all nuclei. It has been proposed elsewhere that the low-energy component of the 
isoscalar dipole mode is determined by surface effects. The existing discrepancies between 
theory and experiment, particularly in the case of the ISGDR in 90 Zr, are significant. The 
resolution of this differences demands substantial effort on both fronts. While extracting 
moments of the distribution will continue to be useful, we suggest that in future studies the 
full distribution of strength be adopted for comparisons between theory and experiment. 

APPENDIX A: THE J w CONTENT OF THE NUCLEAR POLARIZATION 

The first step into the calculation of the RPA response is the computation of the lowest- 
order polarization given in Eq. (^) . Although this step apparently requires the evaluation of 
a six-dimensional integral, the spherical nature of the underlying mean-field potential enables 
one to carry out the four angular integrals analytically leaving a two-dimensional integral to 
be performed numerically. Thus, through a multipole decomposition, the density-dependent 
part of the timelike polarization may be written as [ |36| , |37| 1, 



where all the dynamical information is contained in ITj (g, q'; lo) and the "geometrical" (or 
angular) dependence is given by the function 



oo 




(Al) 



J=0 




(A2) 



M 
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Here D J MX {q) are the Wigner D-functions. Two of these functions may be combined by using 
the following identity [ 36|j3"7 |, 



j dkf&(q, k)P/;,(k, q') = ^L-Sjj.pj^^ q') , (A3) 

so that the three-dimensional integral equation required for the evaluation of the RPA po- 
larization be reduced to a one-dimensional one, albeit one for each angular-momentum 
channel. Computing any specific multipole of the polarization insertion requires the evalua- 
tion of various reduced matrix elements, which are constrained by angular-momentum and 
parity selection rules. Because of the timelike nature of the vertex (7 ) only natural-parity 
states, such as the isoscalar monopole and dipole compressional modes, may be excited. 

Note that there are large computational demands imposed on an RPA calculation of a 
heavy nucleus. As the RPA equations Eq. ( pT|) are solved using standard matrix-inversion 
techniques |j43|, the lowest-order polarization must be computed on every point of a square 



momentum-transfer grid and for every polarization insertion that mixes with U ^. The 
lowest-order polarization must therefore be evaluated several thousands times for a reliable 
extraction of the RPA response. 



APPENDIX: ACKNOWLEDGMENTS 



This work was supported in part by the DOE under Contract No.DE-FG05-92ER40750. 



16 



REFERENCES 



Ronald W. Mayle, Marco Tavani, and James. R. Wilson, Astrophys. J. 418, 398 (1993). 
J.R. Wilson and RW. Mayle, Phys. Rep. 227, 97 (1993). 
S.E. Thorsett and Deepto Chakrabarty, Ast. Jour. 512, 288 (1999). 
Frederick M. Walter and Lynn D. Matthews, Nature 389, 358 (1997). 
P. Danielewicz, Roy A. Lacey, P.-B. Gossiaux, C. Pinkenburg, P. Chung, J. M. Alexan- 
der, and R. L. McGrath, Phys. Rev. Lett. 81, 2438 (1998). 
B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 



C.J. Horowitz and J. Piekarewicz, nucl-th/0010227 



C. J. Horowitz, S.J. Pollock, PA. Souder, and R. Michaels, Phys. Rev. C63, 025501 
(2001); 

R. Michaels, P.A. Souder, and G. Urciuoli, A clean measurement of the neutron skin of 
208 Pb through parity-violating electron scattering, TJNAF experiment E-00-003. 

D. H. Youngblood, CM. Rozsa, J.M. Moss, D.R. Brown, and J.D. Bronson, Phys. Rev. 
Lett. 39, 1188 (1977). 

D. H. Youngblood, P. Bogucki, J. D. Bronson, U. Garg, Y.-W. Lui, and C. M. Rozsa, 
Phys. Rev. C23, 1997 (1981). 

H. P. Morsch, P. Decowski, M. Rogge, P. Turek, L. Zemlo, S.A. Martin, G.P.A. Berg, 
W. Hrlimann, J. Meissburger, and J. G. M. Rmer, Phys. Rev. C28, 1947 (1983). 

G. S. Adams, T.A. Carey, J.B. McClelland, J.M. Moss, S.J. Seestrom-Morris, and 
D.Cook, Phys. Rev. C33, 2054 (1986). 

T.D. Poelhekken, S.K.B. Hesmondhalgh, H.J. Hofmann, A. van der Woude, and M.N. 

Harakeh, Phys. Lett. B278, 423, (1992). 

B.F. Davis et al, Phys. Rev. Lett. 79, 609 (1997). 

H. L. Clark, Y.-W. Lui, D.H. Youngblood, K. Bachtr, U. Garg, M.N. Harakeh, and N. 
Kalantar-Nayestanaki, Nucl. Phys. A649, 57c (1999). 

D.H. Youngblood, H.L. Clark, and Y.-W. Lui, Phys. Rev. Lett. 82, 691 (1999). 

G. Colo, N. Van Gai, P.F. Bortignon, and M.R. Quaglia, [nucl-th/9904051| , 

G. Colo, N. Van Gai, P.F. Bortignon, and M.R. Quaglia, Phys. Lett. B 485, 362 (2000). 

Zhongyu Ma, Nguyen Van Gai, Hiroshi Toki, and Marcelle L'Huillier, Phys. Rev. C55, 

2385 (1997). 

Zhongyu Ma and Nguyen Van Gai, A. Wandelt, D. Vretenar, P. Ring, |nucl-th/9910054 . 

J. Piekarewicz, Phys. Rev. C62, 051304(R) (2000). 

D. Vretenar, A. Wandelt, and P. Ring, Phys. Lett. B 487, 334, (2000). 

J. -P. Blaizot, Nucl. Phys. A649, 61c (1999), and references therein. 



S. Shlomo and A.I. Sanzhur, nucl-th/0011098 



S. Nishizaki, H, Kurasawa, and T. Suzuki, Phys. Lett. B 171, 1, 1986. 
K. Wehrberger and F. Beck, Phys. Rev. C37, 1148, 1988. 

C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 62, 391 (1989); Nucl. Phys. A511, 
461 (1990). 

D. J Thouless, Nucl. Phys. 22, 78 (1961). 

J.F. Dawson and R.J. Furnstahl, Phys. Rev. C42, 2009 (1990). 

A.L. Fetter and J.D. Walecka, "Quantum Theory of Many Particle Systems" (McGraw- 
Hill, New York, 1971). 



17 



[32] B.D. Serot and J.D. Walecka, Adv. in Nucl. Phys. 16, J.W. Negele and E. Vogt, eds. 

(Plenum, N.Y. 1986). 
[33] B.D. Serot and J.D. Walecka, Int. Jour. Mod. Phys. E6, 515 (1997). 
[34] J.D. Walecka, Ann. of Phys. 83, 491 (1974). 

[35] James D. Bjorken and Sidney D. Drell, "Relativistic Quantum Fields" (McGraw-Hill, 
New York, 1965). 

[36] C.J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 62, 391 (1989); Nucl. Phys. A511, 
461 (1990). 

[37] J. Piekarewicz, Nucl. Phys. A511, 487 (1990). 
[38] G.P Lepage, |rmcl-th/9706029i 

[39] Nuclear Physics with Effective Field Theory II, ed. P.F. Bedaque, M.J. Savage, R. Seki, 

and U. van Kolck (World Scientific, 1999). 
[40] R.J. Furnstahl and B.D. Serot, Nucl. Phys. A671, 447 (2000). 
[41] T.N. Buti et al, Phys. Rev. C33, 755 (1986). 

[42] J.R. Shepard, E. Rost, and J.A. McNeil, Phys. Rev. C40, 2320 (1989). 
[43] William H. Press, Saul A. Teukolsky, William T. Vetterling, and Brian. P. Flannery, 
"Numerical Recipes in C" (Cambridge University Press, 1988). 



18 



TABLES 



TABLE I. Various relativistic parameter sets [33]. The scalar mass and k are given in MeV. 



Set 


9 2 s 


9l 


9l 


m s 


K 


A 


L2 


109.63 


190.43 


65.23 


520 








NLB 


94.01 


158.48 


73.00 


510 


800 


10 


NLC 


95.11 


148.93 


74.99 


501 


5000 


-200 



TABLE II. 


Bound single-particle orbitals in 


16 O and low- 


■energy dipole 


(single-particle) tran- 


sitions in three different relativistic 


models. All 


energies are g 


;iven in MeV. 






Orbital 


L2-n 


L2-p 


NLB-n 


NLB-p 


NLC-n 


NLC-p 


is 1 / 2 


41.39 


37.17 


38.75 


34.59 


39.33 


35.18 


lp 3/2 


20.57 


16.68 


19.89 


16.02 


20.77 


16.91 


lp l/2 


12.53 


8.77 


14.10 


10.30 


15.46 


11.65 


l£)5/2 


3.34 




3.44 




4.46 


1.03 


2S 1 / 2 


1.35 




1.55 




2.50 




Transition 


1PV2 _> 25V2 


11.18 


~ 10 


12.50 


~ 11 


12.95 


~ 12 


lp3/2 _^ 1D 5/2 


17.23 


~ 17 


16.45 


~ 16 


16.31 


15.88 


1P3/2 ^ 25 l/2 


19.22 


~ 18 


18.30 


~ 17 


18.27 


~ 17 



TABLE III. Nuclear dependence for the energy of the isoscalar giant monopole resonance in 



three different relativistic models. All energies are given in MeV. 



Model 


16Q 


40 Ca 


90 Zr 


208p b 


L2 


23.2 


27.3 


26.5 


20.1 


NLB 


22.6 


27.9 


24.1 


18.1 


NLC 


21.5 


21.0 


16.9 


13.1 


Exp. 






17.8 ±0.4 


14.2 ± 0.1 
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FIGURES 
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F particles 

FIG. 1. Spectral content of the nucleon propagator in a relativistic Fermi-gas approximation. 




+ 1" 2 + 3" 4 + 5" 6 + 7" 8 

f 



FIG. 2. The real part of (/^IF 10 for the Feynman-Density and Density-Feynman parts of the 
nuclear polarization as a function of the total angular momentum channel. Results are reported 
for 16 and 40 Ca at q = q' = 100 MeV and a; = 10 MeV. In a consistent mean-field formalism these 
quantities should approach ±A, respectively. 
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5 10 15 20 25 

(0 (MeV) 

FIG. 3. The longitudinal polarization for the isoscalar dipole state computed from the timelike 
component of the vector current (solid line) and from the longitudinal component (filled circles). 
In a consistent mean-field formalism — such as the one used here — they should be identical. Note 
that the imaginary component is the longitudinal response. 
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J I I I J 1 — I 1 1 1 I I I I I I I I I I L 



5 10 15 20 

CO (MeV) 

FIG. 4. Isoscalar dipole strength in 16 in lowest-order Hartree (dashed line) and in a con- 
sistent RPA (solid line) approximation. The dot-dashed line is the RPA response with a residual 
interaction that lacks the contribution from the isoscalar component of the photon. The nonlinear 
model NLC was employed in the calculation. 



22 



0.3 



0.25 



i i i r 



i i i r 



i i i r 



i i i r 



Hartree 

RPA 




J I L 



16 0(f =1";T=0) [L2] 
q=300 MeV 
r|=0.05 MeV 




5 10 

co (MeV) 

FIG. 5. Distribution of isoscalar-dipole strength in 16 in a lowest-order Hartree (dashed line) 
and in a consistent RPA (solid line) approximation. The linear model L2 was employed in the 
calculation. 
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0(1 ;T=0) [L2] 
E =7.12 MeV 



E Hal =11.2MeV 
E RPA =8.1 MeV " 




;l I I I | I I I I | I I I I | I I I I | I I I I; 



O(3~;T=0) [L2] 
E =6.13 MeV 



E Har =9.0 MeV 
£^=6.0 MeV 




50 150 250 350 450 550 50 150 250 350 450 550 



q (MeV) q (MeV) 

FIG. 6. Inelastic isoscalar-dipole (left panel) and isoscalar-octupole (right panel) form- factors 
for 16 O in a lowest-order Hartree (dashed line) and in a consistent RPA (solid line) approximation. 
The linear model was employed in the calculation and the experimental data is from Ref. [El]]. 
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FIG. 7. Nuclear dependence of the isoscalar giant monopole resonance in three relativistic 
mean-field models. 
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FIG. 8. Nuclear dependence of the isoscalar giant dipole resonance in three relativistic 
mean-field models. 
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FIG. 9. Low-energy component of the isoscalar dipole strength in three relativistic mean-field 
models. 
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FIG. 10. Momentum-transfer dependence of the low-energy component of the isoscalar dipole 
strength in a lowest-order Hartree (dashed line) and RPA (solid line) approximations. All calcula- 
tions were performed using the NLC set. 
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